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We study effects of interladder coupling on critical magnetic properties of spin ladder systems 
doped with small concentrations of nonmagnetic impurities, using the scaling theory together 
with quantum Monte Carlo (QMC) calculations. Scaling properties in a wide region in the pa- 
rameter space of the impurity concentration x and the interladder coupling are governed by the 
quantum critical point (QCP) of the undoped system for the transition between antiferromagnet- 
ically ordered and spin-gapped phases. This multi-dimensional and strong- coupling region has 
characteristic power-law dependences on x for magnetic properties such as the Neel temperature. 
The relevance of this criticality for understanding experimental results of ladder compounds is 
stressed. 
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Mott insulators on several geometrically designed 
structures of lattices are known to have thcj-spin exci- 
tation gap from the spin singlet ground statetP-. Several 
compounds such as SrCu203B and (VO)2P2O70^ indeed 
have the gap which may be attributed to the underlying 
ladder structure of CU2O3 or V2O3 network. The spin- 
1/2 Heisenberg model on a ladder given by the Hamilto- 
nian 



n 



{i.j) 



o^ • 0-7 



(1) 



with the exchange interaction between nearest-neighbor 
pairs (z, j) describes low-energy spin excitations. The 
spin gap A5 is as large as 0.5 J for the uniform strength 
of J = Ji]a . 

Recent experiments have revealed that the spin ladder 
compound SrCu203 shows surprisingly sensitive depen- 
dence ji^n small concentration of Zn substitution for Cu 
siteg3'B. Even at the doping concentration of Zn as low 
as X = 0.01, the system shows anomalies in the suscepti- 
bility and specific heat at the temperature T ~ 3K. The 
X dependence suggests that the compound changes from 
a spin-gapped insulator at x = to an insulator with an 
antiferromagnetic long-range order (AFLRO) below the 
temperature of the anomaly, T/v, which quickly increases 
with increasing x. In addition, above Tat, it shows a re- 
markable Curie-like behavior x = C /T for the uniform 
susceptibility and a T-linear specific heat Cp. The sensi- 
tive dependence on x and the appearance of the AFLRO 
have a smiilarity to the case of a spin-Peierls compound 
CuGe03Q). 

Since the Zn substitution plays a role to deplete the 
spin- 1/2 local moment on a half- filled Cu dx2_y2 or- 
bital, the ladder model (1) with a dilute depletion of 
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sites has been extensively studied theoretically. Vari- 
ational Monte Carlo, as well as the exact diagonaliza- 
tion results show that a nonmagnetic impurity induces 
the formation of a static spin- 1/2 moment with a sub- 
stantial enhancement of antiferromagnetic correlations 
(AFC) around the impurityB^. Quantum Monte Carlo 
(QMC) results of a singla ladder have supported this re- 
sult of the enhancementEl^ . This enhancement has also 
been suggested in othpi iheoretical approache^for mod- 
els with a spin gapEl'By). The QMC resultJa^ on lad- 
ders have further reproduced the Curie law, x — C jT 
and the linear-T specific heat C^ = 7T with quantita- 
tive agreements of C and 7 with the experimental results, 
when the nonmagnetic impurity concentration x exceeds 
a rather sharp crossover concentration Xc — 0.04-0.05. 
However, below Xc, the local moments induced around 
the depleted sites are so weakly coupled that they do 
not explain the experimental results for 7 and the ap- 
pearance of AFLRO. This is due to the fact that the ef- 
fective coupling between the induced spin- 1/2 moments 
at distance r decreases exponentially as exp[— r/^o] and 
this characteristic energy scale for the typical distance 
between neighboring induced moments becomes far be- 
low the temperature in experiments for realistic choices 
of the AFC length ^0 of a single ladder. On the con- 
trary, the experimental results must be the consequence 
of strongly and antiferromagnetically coupled moments. 
In this letter we study effects of interladder coupling 
to understand this puzzling feature. We show existence 
of four different regions of the critical behavior sepa- 
rated by two crossover lines in the parameter space of 
the interladder coupling Jl and x. Previously employed 
quasi-one-dimensional (QID) and weak-coupling (WC) 
approaches from the single ladder in the literature are 
valid in one of the four regions, whereas we stress that 
a different region among the four is relevant for under- 
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standing experimental indications. This region is gov- 
erned by the universahty of the QCP for undoped sys- 
tems, where the scahng properties are characterized by 
strong couphng (SC) of the induced moments and the 
multi-dimensional (MD) criticality. We discuss the scal- 
ing properties of this antiferromagnetic transition. QMC 
results on coupled ladders as well as experimental indi- 
cations are compared with the scaling predictions. We 
summarize our basic phase diagram for different scal- 
ing properties in Fig.l. Detailed definitions of WC, SC, 
QID and 3D (three-dimensional) or MD are clarified in 
the discussions below. 

First we estimate the scaling properties of the antifer- 
romagnetic transition and the upper bound for the Neel 
temperature in the QID region of weakly coupled (WC) 
ladders. This region is realized when the interladder cou- 
pling is finite but sufficiently weak. In this region, with 
decreasing temperature, ^o within and along the ladder 
becomes much longer than £^l for the interladder direc- 
tion. The effective exchange coupling between two local 
moments induced at two depleted sites each depends ex- 
ponentially on the distance between these depleted sites, 
r, as J{r) o:^ Joe~^^^^ with Jo being the nearest-neighbor 
exchange coupling of the uniform ladder. Here we neglect 
logarithmic corrections. When the impurities are doped 
randomly, the distribution function of the distance be- 
tween two neighboring depleted sites follows the Poisson- 
type distribution Pr{r) = ^e~^^^ for the concentration 
of depleted sites, x = 1/21. Here, the length is scaled by 
the lattice constant of the ladder. From the above rela- 
tions the distribution function for the amplitude of the 
exchange coupling between the moments at neighboring 
depleted sites follows 



P.(J) ^-|(ln|)i4)^. 
l^ Jo J Jo 



(2) 



When T is lowered, only J above the scale of T can con- 
tribute to the growth of AFC. The ratio of such coupling, 
i?, among all the neighboring exchange between depleted 
sites is estimated as 
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where we take Jth ^ T. When x is low enough so that 
/ > ^0 is satisfied, R c^ I - {Jth/Jo)^''^^ follows. There- 
fore the AFC along the ladder is developed over many 
depleted sites only for not too small R and hence only at 



T<Tc, 



Joe" 



(4) 



with a = 7^ 



^. As we have neglected the quan- 



Jo 



tum fluctuation for progressive singlet formation on the 
stronger bonds as is frequently discussed in the random 
exchange antiferromagentB't2P, Tcri based on the above 
argument of the percolation threshold is certainly an up- 
per bound temperature for the growth of a long correla- 
tion length within a ladder. 

When we consider the interladder couplings Jli and 
Jl2, in the other two directions, the necessary condition 
for the growth of 3D correlation over many depleted sites 



in the WC region is similarly given from 



Rc^l-{ 
with 21 = Xjx^l'^ as 
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at small enough x. The upper bound of Tn in the 3D 
region is given by Tcrs- 

From (||) and (^, it turns out that the QID treatment 
is always invalid at small x. The boundary to separate 
the QID and the 3D region is given by 



xd 



VJliJl2 
Jo ' 



(7) 



Below xd^ three dimensional correlations are developed 
before the correlation grows only along the ladder. Equa- 
tions (Q) and (H) show that, in the region of weakly cou- 
pled local moments, the upper bound of Tat is expo- 
nentially small with an essential singular dependence on 
X. In the WC region, the upperbound of the equal-time 
spin structure factor Mirn ^^^ S{Q)/N = ^r{sQS^)e'^^^ 
at the staggered wave vector Q should be scaled as 



N 



lim S{Q)/N = {Msf/3 ex x^ 



(8) 



because only the spin- 1/2 local moments induced around 
the depleted sites contribute to the order parameter. 
Here (Ms) is the staggered magnetization. 

A qualitatively different x dependence of Tn is de- 
rived in the SC region which is governed by the QCP 
at Jli = Jlic and Jl2 = Jl2c- This critical point is 
defined as the interladder coupling beyond which the 
AFLRO is stabilized at T = even without the nonmag- 
netic impurities. To discuss the scaling properties of the 
antiferromagnetic transition in the SC region, we take 
^ = Jlc — JljJl = Jli = Jl2 and Jlc = J Lie = Jl2c 
for simplicity. 

At X = 0, the correlation length exponent v and the 
staggered magnetization exponent [3 are defined as ^ ex 
8~^ and {Ms) ex |^|^ for small \3\. The Neel temperature 
Tn for ^ < is scaled from the critical point as Tn cx 
{—5Y^ where the dynamical exponent z is introduced. 

When X is finite but (5 = 0, we can derive scaling ex- 
ponents for X dependence because x is directly related 
with the inverse cube of the length scale. For example 
Tn should be characterized by the temperature T^^ below 
which the averaged correlation length ^ exceeds the mean 
impurity distance x~^/^. Below T^*^, the local moments 
induced around the depleted sites become strongly cor- 
related antiferromagnetically through the overlap of the 
correlated region of the radius ^. With lowering temper- 
ature, ^ ~ 1/T^ is expected. This gives the scaling 



Tn cx x^ 



(9) 



with p = z/3 in three dimensions. Similarly (Ms) at 
T = is scaled as 

{Ms) cx x^ (10) 

with q = P/du in d spatial dimensions for d > 2. These 
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scaling properties at (5 = govern the SC region. As 
compared to the WC region which shows the essentially 
singular dependence on x, characteristic scales all have 
power-law dependence on x. The crossover boundary 
between the WC and SC regions is given by 

xc = 6" (11) 

with u = dv. Therefore when we fix 5 and increase x, 
physical quantities show a crossover from the scaling in 
the WC region to that in the SC region at this boundary. 

The universality class of the quantum transition at 
X = as a function of J in c^ dimensions is not ex- 
actly known. However, there exists growing support for 
the eqiiivalence to the d + 1-dimensional nonlinear sigma 
modeltSHlS'. The dynamical exponent z is then given by 
z = 1. The universality class of the four-dimensional 
nonlinear sigma model is characterized by the mean field 
exponents z/ = 0.5 and jS = 0.5. Then for the exponents 
in the SC region, we obtain 2p = q = 2/3 and u = 3/2 
for the 3D case. For comparison, we note that, in 2D, 
the long-range order may appear only at T = 0, where 
q ~ 0.26 and u ^ lA arCj-Qbtained because u ^ 0.70 
and f3 c^ 0.36 are expectedE^. It should be noted that 
the boundary Xc may be rather small because u > 1 and 
Xc is extrapolated to around 0.04 even- for the single ID 
ladder as clarified in the QMC resultdS^. 

We next discuss quantitative aspects of the scaling 
properties by using the QMC results. In this letter, we 
show the QMC results of coupled ladders for the con- 
figuration of the two-dimensional square lattice. This 
structure has two basic differences from that of real 
compounds such as SrCu203. One is that the three- 
dimensionality is ignored so that we cannot discuss scal- 
ing properties of T^ directly. The other point is that the 
experimental lattice structure has a certain frustration 
on the ladder plane while it is not considered here. How- 
ever, even with these drawbacks of the present model, it 
is still useful to understand the four different criticalities 
in the WC-QID, WC-MD, SC-QID and SC-MD regions 
with the robustness of the SC-MD region in the sim- 
plest 2D case. More detailed analyses with a 3D realistic 
lattice structure will be discussed elsewhere. QMC cal- 
culations have been performed by using the world line 
algorithm with a clustex^updating methodtSf and a con- 
tinuous time proceduretlP . The system size is N=Lx x Ly 
with Ly being the ladder direction. Several choices of lin- 
ear dimensions are taken for the extrapolation between 
Ly =32 and 128 for Ly/L^ = 4, 6 and 8. 

We first discuss quantitative aspects at x = 0. At 
X = 0, the QCP for the antiferromagnetic transition is 
estimated as Jlc = 0.32 ±0.02 as illustrated in Fig. 2. For 
Jl > Jlci the AFLRO appears at T = and x = 0. Two 
methods have been employed to estimate Jlc- One is to 
measure the Jl dependence of the spin gap for Jl < Jlc 
For finite size lattices, the spin gap A^ is determined from 
the fitting of the uniform susceptibility xo to the form 



Xo(r) = Aexp[-A,/T]. 



(12) 



The limit A^ ^ oo is taken in the scaling of Ag with 
1/N'^. Figure 3 shows the extrapolated values of Ag. 



In the other method, Jl dependence of {Ms) has been 
examined in the region of Jl > Jlc The critical point is 
consistent with Jlc = 0.32 ±0.02 in both of the analyses. 

When X 7^ 0, it is expected that the AFLRO appears 
in all the regions of J^ > 0, if x is not too large. How- 
ever, scaling properties for the x dependence may be 
drastically different for different regions as we mentioned 
above. Here, we show x dependence of (Ms) at Jl = 0.2 
and 0.3. In the following, we show results for regular 
distribution of impurities. The results for the random 
distribution show qualitatively similar behavior although 
the AFLRO is quantitatively suppressed. Detailed anal- 
yses will be discussed elsewhere. Note that at x = 0, 
the gap extrapolated to the infinite size is estimated 
to be As c^ 0.25Jo at Jl = 0.2 while A^ c^ 0.05Jo at 
Jl = 0.3. In comparison with experimental conditions, 
we note thaLihe experimentally observed gap Ag ~ 500K 
in SrCu203E3' with a reasonable estimate Jo — 2000K as 
deduced from the susceptibility of similar but chain- type 
compounds Sr2Cu03 and SrCu02llj^ leads to a similar 
ratio of A5/J0 to the case at Jl = 0.2. At Jl = 0.3, (Ms) 
appears to follow a power-law dependence (Ms) ~ ax^ 
with q ~ 0.26 for smaller x as predicted in the SC-MD re- 
gion of 2D systems. This is reasonable because Jl = 0.3 
is close to the critical point Jlc- The prefactor b in the 
scaling of the crossover boundary defined by Xc ~ bS^-^ 
is roughly estimated as b ~ 0.5 because at Jl = 0.2, 
Xc ~0.02 as we see in 3. Although the numerical results 
are rather preliminary, the results are not inconsistent 
with the scaling prediction. 

We next discuss the relevance of the SC-MD region for 
understanding the experimental results. Because the 3D 
coupling is crucial to realize a nonzero Tn and the 3D 
region dominates over the QID and 2D regions for small 
X as discussed above, the scaling properties in experi- 
mental results for small x need to be analyzed by the 3D 
scaling. Although it is difficult to estimate the exponents 
quantitatively from the available experimental data, at 
least, it is consistent with a quick increase of Tn as a 
function of x with a wide SC region where even x = 0.01 
is involved. Based on this qualitatively promising im- 
plication, further detailed comparisons of experimental 
results and theoretical predictions here would be inter- 
esting. For example, clarification of pressure effects at 
both X = and x 7^ would be desired because rather 
remarkable effects are expected through the control of 
S near the QCP on the basis of the above theoretical 
analysis. In case of CaV205, it has been reported that 
dilute substitutioiupf Ti for V does not make the system 
antiferromagneticE^F . This difference may be understood 
by assuming much smaller interladder exchange coupling 
than in the Cu compound due to presumable orbital or- 
dering of t2g electrons on the V site. 

We also note that the closeness to the QCP has an 
important implication for the mechanism_pf supercon- 
ductivity under doping of itinerant holesElP. Here, we 
discuss two proposals for the mechanism. One is that 
first proposed by one of the authorsE3^ where the carrier 
doping into a spin-gapped Mott insulator may directly 
result in a superconducting state in the region of the 
gap persistence. This proposal was directly applied in 
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the case of ladder modelsclES'. The other is the mecha- 
nism proposed for the high-Tc cuprates, where recovery 
of the suppressed coherence inherent in 2D (or 3D) sys- 
tems nean-ithe Mott transition results in the supercon- 
ductivityEj' E3 . In this mechanism, the closeness to the 
QCP is crucial. The nonmagnetic impurity effects have 
turned out to be a useful probe to measure the distance 
from the QCP of the magnetic order. The experimental 
results appear to support the relevance of the quantum 
critical region and hence support the latter mechanism. 

In summary, effects of interladder coupling on the anti- 
ferromagnetic transition have been investigated for cou- 
pled ladder systems doped with nonmagnetic impurities. 
Critical properties follow different scalings depending on 
four different regions. They are separated by a crossover 
from the weak-coupling to the strong-coupling regions, 
as well as by a crossover from the quasi-one-dimensional 
to the multi-dimensional regions. The strong-coupling 
and three-dimensional scaling behavior governed by the 
quantum critical point of the undoped system dominates 
over a wide region, which leads to characteristic power- 
law dependences of physical quantities such as the Neel 
temperature on x. We propose that the scaling proper- 
ties in this SC-3D region are relevant and important for 
understanding recent experimental indications in ladder 
compounds. 
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Physics, University of Tokyo. 



Fig. 1. Schematic phase diagram for the scahng properties in 
the plane of interladder coupling J^ and nonmagnetic impurity 
concentration x. The quasi-one-dimensional (QID) and multi- 
dimensional (MD) regions are separated by the solid curve x^ 
while the strong coupling (SC) and weak coupling (WC) regions 
are separated by the curve Xc- 



Fig. 2. The spin gap As (triangle) and the staggered magneti- 
zation (Ms) (circle) as a function of Jl scaled by Jo at T = 
and X = 0. Squares and diamonds show (Ms) for T = at 
X = 0.0156 and 0.0208, respectively. The dotted lines are guides 
for the eyes. 
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